On the Coulomb interaction in chiral-invariant one-dimensional electron systems 
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We consider a one-dimensional electron system, suitable for the description of the electronic 
correlations in a metallic carbon nanotube. Renormalization group methods are used to study 
the low-energy behavior of the unscreened Coulomb interaction between currents of well-defined 
chirality. In the limit of a very large number n of subbands we find a strong renormalization of the 
Fermi velocity, reminiscent of a similar phenomenon in the graphite sheet. For small n or sufficiently 
low energy, the Luttinger liquid behavior takes over, with a strong wavefunction renormalization 
leading to a vanishing quasiparticle weight. Our approach is appropriate to study the crossover from 
two-dimensional to one-dimensional behavior in carbon nanotubes of large radius. 
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The recent experimental availability of single-walled fullerene nanotubes has renewed the interest in the study of 
one-dimensional electron systems [[lj-g). In one spatial dimension the Luttinger liquid concept replaces the Fermi 
liquid picture, and provides the paradigm of a system with strong electronic correlations |6| 0. The investigation of 
the generic electronic instabilities of metallic nanotubes has been accomplished in Refs. Jul and There have been 
also recent attempts to look for signatures of Luttinger liquid behavior in single-walled nanotubes that are packed in 
the form of ropes [ ^2[ . Another interesting instance that seems to be feasible from the experimental point of view is 
that of nanotubes in the absence of external screening charges or, at least, with a screening length much larger than 
the typical transverse dimension. The phenomenology of these systems has been studied in Refs. []X3[— p75|] . Anyhow, 
as long as the description of the Luttinger liquid behavior is usually made under the assumption of a short-range 
interaction, the carbon nanotubes with unscreened Coulomb interaction should deserve further theoretical analysis, 
devoted to ascertain possible deviations from the standard Luttinger liquid picture. 

In the particular case of the metallic carbon nanotubes , the Coulomb interaction is also special in that the 
long-range potential does not lead to the hybridization between left- and right-moving electrons in the nanotube. 
At half-filling, the metallic nanotubes have two Fermi points, characterized respectively by the large momenta Kp 
and —Kp, with the typical band structure shown in Fig. [I]. The underlying lattice is such that it allows to arrange 
the modes of the two linear branches at each Fermi point into a Dirac-like spinor, whose components stand for the 
respective amplitudes in the two sublattices of the honeycomb lattice |T^ ] . The kinetic part of the hamiltonian can be 
approximated as H^in = —ta 5k x cti, where a is the lattice spacing and t is the hopping parameter. This means that 
the right-movers have an amplitude ^(r) that alternates sign from one sublattice to the other, while the left- movers 
keep the same sign on both of them. When projecting the Coulomb interaction to the longitudinal dimension of the 
tube, one has to sum over the points of the two transversal rings of the tube C x and C x > 

(9 a {x)Hfp(x') \V(x - x')\ %(x)* s (x')) = *i (r,)*7('i) |— rTT^W)*^^) W 
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It becomes clear that, for a distance x — x' much greater than the lattice spacing, the matrix element is only nonva- 
nishing when a and 7, as well as (3 and 5, have the same (left or right) chirality. 

Thus, the relevant interaction in the study of the metallic carbon nanotubes turns out to be the long-range Coulomb 
interaction between currents of well-defined chirality. In this work we study the genuine effects due to the l/\x\ 
interaction in this kind of one-dimensional systems, focusing on the physics at each of the Fermi points. We disregard 
in this way backscattering processes and residual short-range interactions mixing chiralities, as they have a smaller 
nominal strength (~ 0.1 e 2 /n, in terms of the number n of subbands fl^Jl^j ) and they stay small down to extremely 
low energies (~ t exp(— 75) JT^|). Our present purpose is to discern the low-energy properties of the system and, 
in particular, whether the long-range interaction remains unscreened to arbitrarily large distances. It is known that 
the bosonization approach yields divergent results for some of the observables like, for instance, the plasmon velocity 
at vanishing momentum. This shortcoming is remedied in practice by introducing some infrared cutoff dictated by 
the external environment of the one-dimensional system. Anyhow, it is worthwhile to analyze to what extent the 
singular interaction may be renormalized in the infrared by means of a dynamical screening effect. This problem may 
be also relevant when studying the properties of nanotubes of large radius, since in those systems a crossover from 
two-dimensional to one-dimensional behavior is to be expected when taking the infrared limit. It is known that the 
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Coulomb interaction is strongly renormalized on the graphite sheet [|8[ , what points again at the question about the 
nature of the screening effects in carbon nanotubes with large number of subbands. 

In this Letter we use renormalization group (RG) methods to find the low-energy effective theory of the l/\x\ 
interaction. To be more precise, wc pose the problem of a one-dimensional model with an interaction hamiltonian 
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where ^l(x) and ^r(x) are the electron field operators for the left and the right branch of the linear dispersion 
relation, respectively. The RG method is a sensible approach to deal with this problem since the l/\x\ interaction 
potential (as well as the S(x) potential) gives rise to a marginal four-fermion interaction. The scaling dimension of the 
electron field ^(x) is —1/2 , in length units. This means that the interaction hamiltonian in (||) scales appropriately, 
with a dimcnsionless coupling constant e 2 (in units in which h = c = 1), as the energy scale is reduced down to the 
Fermi level @. 

However, the drawback in dealing with (^) is that it contains a highly nonlocal operator, which makes unclear the 
applicability of RG methods, usually devised to deal with a set of local operators. This problem can be circumvented 
by introducing a local auxiliary field to propagate the Coulomb interaction. The hamiltonian can be written in the 
form 

H = iv F J dx (<S>+(x)d x y R (x) - V+{x)d x V L (x)) +ejdx (V+(x)* L {x) + <S>+(x)Vr(x)) <f>(x) (3) 

where the 4>{x) field propagates the interaction 

i{T<f>{x, t) c/>(x',t')) = ±6(t - f)j^- . (4) 

We may think of 4>(x) as the scalar potential in three-dimensional quantum electrodynamics. However, the differ- 
ences with that theory in the present case are notorious since the propagation of (f>(x) is that of a genuine field in three 
spatial dimensions, while the electrons are confined to one dimension. In general, one may expect a better infrared 
behavior in the present model. The propagator of the 4>{x) field in momentum space can be read from the relativistic 
expression p0| , after sending the speed of light to infinity, 

,<T*(*. t) *{x ,t)) = J ^ J w q2 + q2y + q ,_ le (5) 

The usual one-dimensional propagator ~ log(|g|/A) is recovered from (|^) upon integration of the dummy variables 
q y and q z . We remark that the ultraviolet cutoff A for excitations along the y and z transverse directions is needed 
when projecting the three-dimensional interaction down to the one-dimensional system. 

The usefulness of the representation (||) can be appreciated in the renormalization of the model at the one-loop 
level. We study the scaling behavior of the irreducible functions as the bandwidth cutoff E c is sent towards the Fermi 
level, E c — > 0. The self-energy to the one-loop order is 
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The limit k — > has to be taken carefully in this expression, by first combining the two denominators with the use of 
Feynman parameters (2(| . Finally we get 

i e 2 f 1 1 f E " dp f°° dujp vpk 



i S(fc,0) = — — / du^ / £ 
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The term linear in k in E(fc,0) represents a renormalization of the Fermi velocity, which grows upon integration of 
the high-energy modes. Obviously, there is no correction linear in u>k renormalizing the electron wavefunction at the 
one-loop level. This is consistent with the fact that the integration of high-energy modes at ~ E c does not renormalize 
the three-point vertex T . We stress the difference of the logarithmic renormalization of Vp with respect to the usual 
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finite corrections due to a short-range interaction. The nontrivial scaling of vp is a genuine effect of the long-range 
Coulomb interaction, which also takes place in higher dimensions [ [l8|]2l| l. 

Incidentally, the above computation exemplifies how the Ward identity that ensures the integrability of the Luttinger 
model does not hold in the present case. The Ward identity is a relation between the electron Green function G(p, u> p ) 
and the three-point vertex T(p, ui p ; k,u)k) at a given branch For the right-handed modes, for instance, it is 

T{p, uj p ; k, Lo k ) = £ : £ (8) 

oJk~ vpk 

By focusing on the singular dependences on the bandwidth cutoff E c , one can check that (||) is already violated in 
our model to first order in perturbation theory. Actually, the dependence of the vertex T on the variables {k,ojk) of 
the external interaction line is 

„, , . e 2 k log A 
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We notice that the Ward identity would be satisfied if the scaling could be implemented simultaneously in the 
transverse ultraviolet cutoff A and the bandwidth cutoff E c , that is, by taking A = E c . However, in a real system 
the scaling in A gets locked by the finite cross section of the wire, while only the scaling in the longitudinal direction 
operated by E c is allowed. In this respect, the gauge invariance of quantum electrodynamics is broken by the anisotropy 
of the electron system, as felt by the propagation of the three-dimensional electromagnetic field. 

The renormalization of vp at the one-loop level is not, in general, a sensible effect from the physical point of view, 
since the propagator of the 4>(x) field is drastically modified by the quantum corrections. In what follows we implement 
a GW approximation in order to take into account the dynamical screening due to plasmons. The suitability of this 
approximation in the study of one-dimensional systems has been recently shown in Ref . . The same approach has 
been also tested in the study of the crossover from Fermi liquid to Luttinger liquid behavior |22j ] , as well as in the 
study of singular interactions in dimension I < d < 2 |24|,^5| . 

The self-energy U(k,ujk) of the <j)(x) field is given at the one-loop level by the sum of particle-hole diagrams with 
modes in the left and the right branch of the dispersion relation. The number of contributions depends on the number 
n of different Dirac fermions, so that 

m(fc,^)=m- . |fc2 _^ (10) 

In the case of carbon nanotubes we have n = 4, taking into account the two Fermi points and the spin degeneracy 
, but it is also conceivable that in the process of renormalization n is given effectively by the number of subbands 
within the energy cutoff E c . In the spirit of the GW approximation, we consider vp as a free parameter that has to 
match the Fermi velocity in the fermion propagator after self-energy corrections. We recall that the result ( |l0| ) turns 
out to be the exact polarization operator in the model with short-range interactions, with an unrenormalized Fermi 
velocity v F — v F In the present case, though, the violation of the Ward identity (||) already signals that we do 
not have at work the precise cancellation between electron self-energy insertions and vertex corrections characteristic 
of the Luttinger model. The relevant remnant is actually the renormalization (Q) of vp , and it can be taken into 
account self-consistently by replacing vp in (|l^) by the renormalized value of the Fermi velocity. 
The particle-hole processes lead to a modified propagator of the 4>{x) field 

mk,u) = i/ (- log( |^ /A) + n ( fc ^)) ( n ) 

The expression ( |T]] ) provides a sensible approximation for the scalar propagator, as it incorporates the effect of 
plasmons in the model. Thus, our approach is that of using the scalar propagator ( |TT] ) in the renormalization of 
the Fermi velocity and the electron wavcfunction. We compute the electron self-energy by replacing the Coulomb 
potential by the dressed interaction ( |fl| ) 

27r J- Ec ^ J-oo 2tt i(wp + Wfc ) - v F (p + k)\- n ^ 1 ^L_ log(M/A) 

Alternatively, one can think of the diagrammatics encoded in Eq. (|T|) as the leading order in a 1/n expansion, in 
a model with n different electron flavors. It can be shown that this approximation reproduces the exact anomalous 
dimension of the electron field in the Luttinger model with conventional short-range interaction p7[. In our case, 
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such approximation to the self-energy is also justified since it takes into account, at each level in perturbation theory, 
the most singular contribution at small momentum transfer of the interaction. Due to the cancellation of fermion 
loops with more than two interaction vertices which still takes place in the same way as in the Luttinger model, the 
representation ( |l2|) for the self-energy only misses the effects of vertex and electron self-energy corrections, which may 
be incorporated consistently in the RG framework by an appropriate scaling of the renormalized parameters. 
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The only contributions in (12) depending on the bandwidth cutoff are terms linear in u>k and k. In this respect, 
it is worth mentioning that, although the usual perturbative approach gives rise to poles of the form k 2 /(u>k — vfR) 
in the self-energy Pq| , these do not arise in the GW approximation. There is no infrared catastrophe at u>k ~ vpk, 
because of the correction in the slope of the plasmon dispersion relation with respect to its bare value vf- The result 
that we get for the renormalized electron propagator is 

G _1 (/s, uJk) = Z^, 1 (u k - v F k) - £(fc,cjfc) 
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where f(p) = 1 — ne 2 log(|p|/A)/(27r 2 WF) , r = vf/vf and z\j 2 is the scale of the bare electron field compared to 
that of the cutoff-independent electron field 

Vbare(Ec) = Z^** . (14) 

In the RG approach, we require the cutoff-independence of the renormalized Green function, since this object leads 
to observable quantities in the quantum theory. For this purpose, the quantities Z-% and vp have to be promoted to 
cutoff-dependent effective parameters, that reflect the behavior of the quantum theory as E c — > and more states are 
integrated out from high-energy shells of the band. Regarding the problem of self-consistency for the renormalized 
value tip of the Fermi velocity, we find two possible solutions leading to different physical pictures: 

i) large-n limit solution. In this limit we know that the expression (|Io| ) gives the exact result for the polarization 
operator, since any vertex or self-energy correction makes any diagram subdominant from the point of view of the 
1/n expansion. Then the correct choice for vp has to be the fixed-point value of the Fermi velocity. Self-consistency 
is therefore attained by requiring that the solution of the scaling equation 
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matches the fixed- point vp in the limit E c — > 0. It can be checked, however, that any finite value of vf in the 
right-hand-side of ( |l5| ) does not lead to a strong enough flow to reach vp at E c = 0. The only self-consistent solution 
is found for vp = oo. For this value the right-hand-side of ( |l5| ) has a finite limit, which produces the asymptotic 
scaling of the Fermi velocity Vp(E c ) ~ — e 2 /(127r 2 ) \og(E c ). 

ii) one-band solution. If we stick to the picture in which we only pay attention to the left and right linear branches 
of the dispersion relation, we have to assume that (|l^) only provides an approximate expression for the polarization 
operator. The incomplete cancellation between electron self-energy and vertex corrections to that object arises from 
the mismatch between E c and the transverse cutoff A. As a consequence of that, the renormalized Fermi velocity 
in ( |Io| ) gets an effective dependence on the variable log(£' c /A) and vp is to be taken as the scale dependent Fermi 
velocity, vp = vp(E c ). With regard to the carbon nanotubes, this is consistent with the regime in which the scaling 
has progressed to distances larger than the radius of the nanotube. The RG flow equations then turn out to be 
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E c — log Z*{E C ) = v - - y (16) 
dE c Sy/f(E c ) 
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Eq. ( p"7| ) is now the requirement of self-consistency, with vp = vf(E c ). As mentioned before, the three-point vertex 
only gets the cutoff dependence given by the wavefunction renormalization in (E6h. This means that the electron 
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charge is not renormalized at low energies in our field theory framework. The behavior of the effective interaction is 
therefore completely encoded in Eq. (fl7|), which can be rewritten for the effective coupling constant g = e 2 /(4tt 2 vf) 
in the form 

E °w. 9iEc) = I + mk^) (18) 



The right-hand-side of Eq. (|18|) vanishes as e /(An vf) — * and it could still admit a solution of the form g 



—go/log E c , but this is not realized in the present regime as the equation 850 = \/l + go — 4/3 + 1/(3-*/ (f + 50) 3 ) 
does not have any real solution. The flow of g(E c ) given by ( |l8| ) quickly approaches some fixed-point value, after 
which it becomes little sensitive to further scaling in the infrared p8[ . A plot of the flow for different values of the 
bare coupling constant at E c = A is given in Fig. @, which includes for comparison the scaling behavior of the large-n 
solution. 

The existence of two different regimes in the model is not surprising, since the dependence on the number n of 
subbands is the way in which the system keeps memory of the finite transverse dimension in the carbon nanotube. A 
nanotube of very large radius, for instance, leads to a picture in which a large number of subbands are stacked above 
and below the linear branches in Fig. [|. This system falls into the description of the model with a very large n value. 
On the other hand, there is always a sufficiently small value of E c for which all the subbands, but those contributing 
to the gapless part of the spectrum, become higher in energy than the RG cutoff. At that point, one proceeds paying 
attention to the linear branches crossing at the Fermi points alone, what leads to the one-band regime of the model. 

^From the point of view of the real space, the above transition in the number n of subbands represents the crossover 
from the two-dimensional regime to the effective one-dimensional description of a nanotube of large radius. Actually, 
in the limit of a very large n we recover the scaling behavior of the effective coupling constant in the graphite layer, 
which flows towards the free fixed-point. This scaling of the coupling constant is actually what makes consistent the 
linear quasiparticle decay rate recently measured j2{| with the metallic properties of graphite |30];[3l| . In a carbon 
nanotube, the scaling appropriate for the graphite regime has to be followed up to a distance scale of the order of 
the diameter of the nanotube. The value of the coupling constant at that scale is what dictates the bare value for 
the one-dimensional regime ii). The one-dimensional effective field theory contains the explicit dependence on the 
ultraviolet cutoff A, that is actually needed to fix the conditions at the crossover between the two regimes. 

The issue of the renormalization of vf we have discussed is important, since it implies a reduction in the strength 
of additional interactions in the system, whether short or long-range, as the effective couplings are all given by the 
couplings in the interaction hamiltonian divided by vp ■ We stress that this renormalization of the Fermi velocity 
is the relevant effect for a phenomenology at realistic energies, because the short-range interactions like umklapp or 
backscattering are nominally tiny and their RG flow only becomes appreciable at extremely small energies. 

Incidentally, such a renormalization of vp affecting every effective interaction in the model may also be present in 
small chains, contributing to explain the insulator-metal transition by the effect of the Coulomb interaction observed 
in the exact diagonalization of finite rings |3^ ]. The study of low-dimensional systems carried out in Ref. |33| also 
indicates a reduction of the electron correlations similar to our findings in the RG framework. 

The reason for the distinctive behavior of the Coulomb interaction in the fullerene tubules is the preservation 
of the chiral invariance in such systems. This comes from the different symmetry properties that characterize the 
left and the right modes at each Fermi point, leading to the selection rule that forbids their hybridization by the 
interaction. Thus, although it has been proposed that in a generic one-dimensional electron system the 2kp and 



higher harmonics may lead to strong correlation effects in the presence of the Coulomb interaction 34 1 , the absence 
of such hybridization makes the carbon nanotubes fall into the more conventional Luttinger liquid regime. As stated 
above, the wavefunction renormalization ( fl6|) is just the generalization for the long-range interaction of the expression 
giving the anomalous dimension of the electron field in the Luttinger model. The important point is the existence 
of a stable fixed-point for the effective coupling constant in the infrared. All the properties of the Luttinger liquid 
universality class are then expected to hold in the metallic carbon nanotubes, with a sensible renormalization of the 
effective interaction strength depending on the radius of the tubule. 
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FIG. 2. Flow of the effective coupling constant for different bare values, computed in the large-n solution (thick lines) and 
in the model with only one band (thin lines). 



8 



